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On a Method of Reducing Occultations of Stars by the Moon; 
together with the Reduction of Occultations observed on Three 
Occasions at the Liverpool Observatory. By H. C. Plummer, 
M.A. 

The present paper is the result of an undertaking to reduce 
three series of observations of occultations and to render them 
available for the discussion of the Moon’s semi-diameter. As I 
was not acquainted with a method which did not seem to leave 
something to be desired in point of clearness or convenience, I 
adopted one which offers certain considerable advantages. This 
method seems to differ essentially from those which have been 
previously employed, and it therefore appears desirable to 
describe it in some detail. The first object of the theory is to 
obtain a quite general equation expressing the geometrical con¬ 
ditions, on the assumption that the figure of the Moon is a sphere 
and that the line joining the observer to the star is a tangent 
line to the Moon at the time of the observation. 

Let {x, y, z) be the position of the observer ; (X, Y, Z) of the 
centre of the Moon ; and ( l , m, n) direction cosines of the star, 
referred to rectangular geocentric axes. Let R be the radius of 
the Moon. Then 

i—x _ n—y _ Z—z 

l m n 

l—x _ v—y Z—z 

X —x Y— y Z —z 

represent the lines joining the observer to the star and to the 
centre of the Moon respectively. If, now, 'J is the angle between 
these lines, 

l(X-x) + m(Y-y) + n(Z-z) 

=cos 6 . {(X-x ) 2 + (Y-y ) 2 + (Z-z) 2 } ; 

R=sin-J. s/{(X-z) 2 + (Y -y)z + (Z-zf) ; 

(Z(X— x) + m(Y—y) -f n{Z—z)} 2 + R 2 

=(X-xy + (Y-yy + (Z-z) 2 . 

Xow let p be the Earth’s radius and <]> the geocentric latitude 
at the place of observation ; T, the sidereal time (or R.A. of the 
zenith); A, D, the apparent R.A. and declination of the star; 
«, c, the geocentric R.A. and declination of the Moon. Also let 
r, P, s be the Moon’s radius vector, parallax, and semi-diameter 
respectively. Then, if the axis of x has been taken through the 
Eirst Point of Aries and the axis of z through the Pole, 

X=rcos8cosa; Y=rcosSsina; Z=rsin8; 
£=cosDcosA; m=cosDsinA; w=sinD; 
x=p cos 0 cos T; y=jocos^sinT ; z=|osin0. 
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Making these substitutions in the above equation, we obtain 

R, 2 + (r [cos 8 cos D cos (A — a) + sin D sin 2 ] 

—p[cos D cos cos(A — T) + sin D sin <^]} 2 
= r 2 +p 2 —2 rp (cos <5 cos <j> cos (a — T) + sin 8 sin . 

Now 

cos § cos D cos (A—«) + sin D sin 8=:cos A ... (1) 

cos D cos ^ cos (A—T)H-sinDsin0==cosZ ... (2) 

cos 8 cos <p cos (a —Tj + sin 8 sin— <f>—cos Z' ... (3) 

where A is the geocentric distance of the star from the Moon’s 
centre, Z is the true Z.D. of the star and Z ; is the true Z.D. of 
the Moon’s centre. Hence 

R 2 + (r cos A— p cos Z) 2 =r 2 4 -p 2 — 2 rp cos Z' 
or R 2 =r 2 sin 2 A + p 2 sin 2 Z—2rp(cos Z' — cos A cos Z). 

But r sin P— 1 ; rsins=R; the Earth’s equatorial radius being 
unity. 

Therefore 

sin 2 s=sin 2 A+p 2 sin 2 P sin 2 Z—2p sin P (cos Z '—cos A cos Z) (4) 

and this is the geometrical equation sought. 

We are now in a position to derive the equations of condition. 
Since p and. P do not occur independently in the fundamental 
equation, p may be written for p sin P. Also put S equal to 
sin 2 s, and let 

C=cos Z' —cos A cos Z ... ... ... (5) 

2 =sin 2 A+p 2 sin 2 Z— 2pC ... ... (6) 

Then, if the observed and tabular quantities involved were quite 
accurate, the equation 2 would be satisfied exactly. Actually, 
however, the quantities s, p, a, 8, A, D, and t, the mean time of 
the observation, require corrections ds, dp, da, do, dA, dD, and 
dt. Substituting s + for s, &c., an equation is obtained of the 
form : 

2—S=sin 2 s . ds—'Bdp—Rda — P/dA — Fc?o — FWD —Hdt (7) 


wherein 


02 

E= — • F= 

da’ 


?? . . uv_ . 

08 ’ 0 A ’ 0 D * 


H= ^2 = 02 tfT + 02 
dt 0 T' dt ' 0 a 


da 02 dS 
dt 08 dt 


=G-+E~ a - +F — 
dt dt 


... ( 8 ) 
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and 

B=^ = 2 (psin3Z-C). 

op 

Hence 


147 


... (9) 


E=—2 cos A (cos A )—‘Ip i UicosZ')—cosZ —(cos A )1 

0 a l ca 0a i 

= — 2 {(cos A— p cos Z) cos c cos D sin (A—a) 

— p cos c cos ^ sin (a—T)} ...(10) 

0 0 
F= —2 cos A — (cos A) — 2p < --(cos Z)—cosZ —(cos A) V 
do l cc do J 

——2 {(cos'A —p cos Z)[ —sin c cos D cos (A—a)-}-sinD cos 5 ] 
+p[—sine cos ^ cos (a — T) + cos £sin^>]} ... (n) 


and 


02 dT dT r 
G= ST • *=* ' {-^cosZ, T (cosZ) 


— 2 p 


0 0 

5l ,(cos Z )-cos Ag^(cosZ)J J- 


dT 


= — 2 'p—j- {{p cos Z—cos A) cos D cos <p sin (A—T) 

+ cos c cos (f> sin (a—T)} ... (12) 

Finally, 

0 f 0 

E'= — 2 cosA^— (cos A) + 2 p { cos Arr— (cos Z) 

O-ix L C/-i-X 

0 1 

+ cos Z g^( cos r 

= 2 {(cos A— p cos Z) cos c cos D sin (A— a) 

— p cos A cos D cos f sin (A—T)} ...(13) 


and 


a r a 

F'=—2 cos A — (cos A) + 2p <j cosA~ (cos Z) 


3 1 

+cosZ —(cosA) j 


=2 {(cos A —p cos Z)[cos l sin D cos (A—ct)—cos D sin S] 

-f-^ cos A[—sin D cos <j> cos (A—T) + cos D sin (f]\ ...(14) 


We thus have the means of calculating accurately all the co¬ 
efficients of the equation of condition. It is only necessary to 
add, since p is not itself a tabular quantity, that 

dp—p cos P. c£P+sin P . dp ... ... (15) 
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It is now necessary to reconsider the above formulae from the 
point of view of the computer. In the first place cos Z and 
cos Z' can be found by means of (2) and (3). It is assumed 
that 7-figure logarithms are used, and it may be pointed out that 
perfect accuracy will be attained if the terms of 2, none of which 
in fact exceeds '0005, are calculated to ten places of decimals. 
Equation (1) is unsuitable for obtaining sin 2 A, but from it can 
be derived 

sin 2 A=i—cos 2 A = i— {cos(D—8) — 2 cosS cos D sin 2 |-(A —o)} 2 

=sin 2 (c —D)+4 cos 8 cos .D cos (S -D) sin 2 1- (A—a) 

—4 cos 2 S cos 2 D sin 4 -§■ (A—a) 

=sin 2 (8—D) + cos D cos c cos (S— D) sin 2 (A —n) 

-f sin 2D sin 28 sin-* * (A—a) ... ... ... (16) 

In this form the last term is generally negligible, but in any 
case it is easily added as a small correction; its importance 
depends only on the fact that it is always positive. From 
sin 2 A can be deduced cos A, and then we have C from (5). This 
is the most convenient way of calculating 0, but it must be said 
that the resulting value of 2pC cannot be trustworthy beyond 
nine places, since 0 is known only to seven places of decimals. 
This falls short of the highest standard of accuracy, but I have 
generally contented myself with this method : partly because of 
its convenience, partly because the error introduced is an 
accidental one, small in comparison with those which depend on 
the star’s place and the observed time, and partly because other 
methods failed to reveal errors in the simpler process, carefully 
performed, such as seemed to necessitate or warrant the extra 
labour involved by their adoption. Two alternatives may, 
however, be mentioned here. The angles Z and Z' are known 
from their cosines. Then (5) may be written 

C=2 sin ■§-(Z-f-Z') sin i (Z —Z')+2 cos Z sin 2 ■§-A ... (17) 

Or referring to the geocentric spherical triangle formed by the 
zenith (Z), the Moon’s centre (M) and the star (S), it is seen 
that 

C=cos Z r —cos A cos Z=sin A sin Z cos ZSM. 

But if P is the pole, ZSM=PSM+ZSP='sr-f say. 

C=sin A sin Z cos (w+£) ] 

where sin w=cos 8 sin (A—a) / sin A i... ... (18) 

sin 4=cos <f> sin (T—A) / sin Z j 

These latter equations may not give the auxiliary angles without 
ambiguity : in that case reference to a figure will remove the 
difficulty. The use of equations (18) involves considerable extra 
labour, but they possess the advantage that a comparison of the 
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value of C so obtained with the value obtained directly from (5) 
furnishes a very complete if not a very accurate check on the 
work. There is no further difficulty in finding S — S, which is at 
present expressed in circular measure. 

In finding the coefficients of equation (7) 5-figure logarithms 
are of course amply sufficient. In my own computations, B has 
been found from equation (9). In finding E, F and G the 
following slightly modified approximate for mu lie have been 
used :— 

E= — 2 cos S {(1 —p cos Z) cos D sin (A —a) 

—p cos (f» sin (a —T)} ... (19) 

F= — 2 {(1 —p cos Z) sin (D — c ) + ^[sin <f> cos 3 

—sine cos $ cos (T—a)]} ... (20) 

dT 

G= 2 p cos 0 cos D {sin (A —a) cos (A—T) 

—p cos Z sin (A — T)} ... (21) 


cl clS 

In equation (8) — and — are most easily expressed in seconds 
d/t do 


of arc per second of mean time; for the sake of uniformity 
therefore we take log ^2 = 1 ’47831. It is easy to see, either 


by comparing equations (10) and (13), and (11) and (14), or from 
general considerations, that E'=—E and F ; = — F approxi¬ 
mately. In consequence of this it has not been thought worth 
while to calculate the coefficients E' and F' separately. In fact 
these coefficients are not needed with great accuracy, for if con¬ 
siderable values be assigned to dA. and d\), terms of the second 
order will become necessary in equation (7). The term in dp is 
most convenient if we are dealing with observations made at a 
single observatory. If we have to combine observations made at 
different observatories, it will be necessary to remove the term 
in dp in equation (15) and to use —Bp cos P (which differs little 
from — B) as the coefficient of a term in dP on the right-hand 
side of equation (7). If it is desired to improve p by adding dp, the 
quantity B sin P dp must be added to the left-hand side of the same 
equation. For an increase of altitude of 209 feet, dp=' 00001, 
while B sin P is not likely to exceed - oooi. In the final equation 
we divide 2 — S by sin 1" in order to transform the residual 
from circular measure into seconds of arc, and divide the whole 
equation by sin 2 s so as to make the coefficient of ds unity. 

It is well known that refraction exercises an effect, which is 
of a parallactic nature, on the circumstances of an occultation. 
If there were no terrestrial atmosphere it would be necessary for 
the observer to change his position in order to see the occultation 
at the same time that he makes his observation under actual 
conditions. The effect, which is generally evanescent and always 
small, is presumably neglected in the Greenwich reductions, since 
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do reference is made to it, but it is taken into account by those 
who follow Bessel, as for instance by Kiistner in his admirable 
memoir on “ Nine Occultations of the Pleiades ” It is easy to see 
that the effect is allowed for, if the radius p is increased by an 
amount x such that 

1 +x=fi sin 2/sin Z 

where y is the refractive index of the air at the surface of the 
Earth, and 2 is the apparent Z.D. of the star (cf. Chauvenet’s 
Astronomy , 5th ed., voi. i., p. 516). Now z—Z — a tan Z where 
e tan Z is small. Consequently 

i+aj=yu(sin Z — a sin Z)/sin Z=p(i — a) 
car x=y —1— ya—y—i — a. 

It is quite clear that an accurate calculation of this quantity 
can be made only if y. and a are known with great exactness. 
If we adopt Gladstone and Dale’s law as to the dependence of 
the refractive index on pressure, and Mascart’s numbers (Comptes 
Mendus, t. lxxviii., p. 617) for the index itself and the temperature 
coefficient, we have 

y —i = *ooo2923(i — *00383^)— . 

P° 

And using Bessel’s formula and numbers for the refraction we 
have 

a=a/3 A y A =a ^ I+Ml (1— et) 1+n 

where m and n are small, and «=‘00364. At the zenith, when t— 
©°C and 77=760 mm., a=*ooo2927, we may put then <*=‘0002927 
—a for any Z.D., and let pjp Q =i-\-q. If we make these trans¬ 
formations in the expression for x, we obtain a formula which is 
certainly erroneous in a slight degree, inasmuch as it gives a 
negative value of x at the zenith when q—t=o. This may be 
due to the fact that Bessel’s temperature coefficient ‘00364 is too 
small. However it is only necessary to retain terms amounting 
to *000001, and from an examination in detail it is concluded 
that we may take x—a\i +<7—‘00364^), and that no serious error 
will be introduced by assuming, as has been generally done, a 
mean state of the atmosphere, in which case we find the fol¬ 
lowing :— 

Z= o° 50° 6o° 65° 70° 75 0 8o° 85° 

a?=*oooooo 000 001 002 003 005 012 038 

It is to be feared that the constants of refraction may be in¬ 
sufficiently exact in the only case where they are of importance 
here, i.e. at great zenith distances. For this reason the correction 
has not been applied in the reductions which follow, but its 
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magnitude in the few instances in which it amounts to o" oi may 
be recorded here, though, as a matter of fact, in no case has 
it any importance. In the final equations the coefficient of 
dp= —B/sin 2S=B' say, and the residual is (S, —S)/sin 2s sin i n . 
It is only necessary then to add B sin P . xj sin 2 s sin 1"= —B'Prc, 
where P is expressed in seconds, to correct the final equation for 
the effect of refraction. The equations affected are confined to 
Series I, and the corrections are — o"'02 in No. 3, —o"'oi in 
Nos. 5, 6, 7, and 9, and +o //- oi in No. 11. 

The accompanying reductions refer to observations made by 
the Director of the Liverpool Observatory on the following 
dates :— 

1897 July 23 (Monthly Notices, vol. lviii., p. 20) 

1898 Dec. 27 ( „ „ lix., p. 165) 

1899 Dec. 16 ( „ „ lx., p. 222) 

On the first of these dates the Pleiades were occulted, and on the 
other two the Moon was eclipsed. A few corrections have been 
made in the designations of the stars, and in several instances it 
has not been found possible to identify the stars occulted. In 
Series I the mean places of the stars have been taken from 
Elkin’s memoir on the Pleiades, except in those cases where the 
apparent places were given in the Nautical Almanac of the year. 
In columns 5 and 6 are to be found apparent coordinates (A andD) 
as used in the reduction. The stars in Series II are all taken 
directly from a list prepared by the photographic method at the 
Poulkova Observatory (Astr. Nachr., No. 3533). For Series III 
the three stars B.D. + 22 0 , Nos. 999, 1004 and 1006, were taken 
from the Astr. Ges. Catalogue (Berlin). The coordinates of 
No. 1011 are the result of my own reduction of a plate taken at 
the University Observatory, Oxford, by Mr. Bellamy, to whom, 
as well as to Professor Turner, I am indebted for kind assistance 
in this connection. The corrections deduced from the plate to 
the Berlin meridian places are these : — 

o s u 

B.D.+ 22, 999 AA=—C05 AD=—o*4 

IOO4 -f- 0'02 —0*2 

1006 +0'02 +0'2 


As regards the places of the Moon (a and t), columns 7 and 8 
contain the values deduced by interpolation from the Nautical 
Almanac with the following assumed corrections added :— 


1897 July 23 

1898 Dec. 27 

1899 ® ec - *6 


S 

Aa= +0*06 
+ 0*20 
+0-38 


Ac= — 2 ‘o 
+ 2*3 
-{- o*8 
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These are based on observations made at Greenwich, in the 
third case directly, and in the two former by estimation, not 
very reliable probably, from the results of previous and subsequent 
nights. The unpublished results were kindly communicated to 
me by the Astronomer Royal. The assumed position of the 
observatory is that given in the Nautical Almanac , and no cor¬ 
rection has been made for the altitude, which is about *000011 in 
units of the Earth’s radius. The omission contributes to the 
term dp, and might be allowed for in the way already explained. 
But this is unnecessary if the term is eliminated from the final 
equations, and if this is done the result will be independent of 
the assumed figure of the Earth. 
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I. 

Occultation of the Pleiades , 1897 July 23. 


i. 

+ 1*2 5 = ds 

— •0241 dp 

+ 7364 da 

+ -5833 d8 

+ '6497 dt 

2. 

H-4‘l7 = ds 

+ ’0248 dp 

+ 7511 da 

+ ’5680 dS 

+ '6403 dt 

3- 

+ 2‘02 = ds 

+ ‘8683 dp 

4* ‘8256 da 

-•4379 d5 

+ •5811 dt 

4- 

+ 206 = dt 

-•1375 dp 

+ *6379 da 

+ 7121 dS 

+ •5861 dt 

5- 

+ 7-11 = ds 

+ *5206 dp 

+ *9125 da 

+ •1109 dS 

+ -6690 dt 

6. 

- 3*33 = ^ s 

+ *5895 dp 

4- *9047 da 

— •0007 dS 

+ -6320 dt 

7- 

— 2-03-ds 

4- *6566 dp 

4- '9065 da 

-•0739 ds 

+ -6439 dt 

8. 

-r 5*14= ds 

4- *1032 dp 

4-7607 c/a 

+ -5634 dS 

+ -6043 dt 

9- 

+ 0*45 = ds 

4- *8544 dp 

4- *8055 da 

— •4760 ds 

+ •3115 dt 

io. 

4- 3-65 = ds 

-r *2025 dp 

4- *8065 da 

+ *4623 dS 

+ -5996 dt 

ii. 

-5*08 = ds 

— ■5788 dp 

— ‘904I da 

— •0263 dS 

-•3391 dt 

12. 

— 4*40 = c/s 

— *0049 dp 

— *6912 </a 

— •6449 <78 

— •2701 dt 


II. Occultations during the Lunar Eclipse , 

1898 December 27, 

I. 

4- 5*63 = £?s 

— '4064 dp 

+ 7122 da 

+ •9433 d5 

+ ‘0961 dt 

2. 

+ 2'85 = ds 

+ •1580 dp 

+ *9140 da 

+ •1099 dS 

+ •3881 dt 

3- 

4- 2'20 = ds 

+ •1345 dp 

4* 9130 da 

+ •1113 dS 

+ •3814 dt 

4- 

+ 4’9S = ds 

+ -5078 dp 

+ *4789 cfa 

— •8603 dS 

+-2416 dt 

5- 

—o - 8o = ds 

-•4159 dp 

+ '3996 da 

+-8985 dS 

+ *1036 dt 

6. 

+ 080=ds 

+ -2968 dp 

4*‘8192 c/a 

-•4530 d8 

+ •3574 dt 

7- 

+ i-64 = c?s 

— •0129 dp 

4-*9181 da 

+ ‘0419 dS 

4- *3480 dt 

8. 

+ 5-96 = ds 

+ •0157 dp 

4- *9232 c/a 

— •0189 dS 

4- *3538 dt 

9- 

—121 =ds 

+ •2531 dp 

— •6465 da 

— 7074 dS 

— •4214 dt 

10. 

—014 -ds 

—3473 # 

— *8073 c/a 

+ '473 s ^8 

— *3222 dt 

ii. 

— 278 =ds 

-‘4797 

- *0493 da 

-r ’9946 d5 

— *0810 dt 

12. 

— 232 = ds 

— *2240 

— *8311 da 

+ ‘4265 dS 

-•3449 dt 

13- 

— 052 = ds 

+ •3234 dp 

— 7184 c/a 

— •6216 d5 

— •2316 dt 

14. 

+ 3-43 = « s 

— *2108 

— ’8273 da 

+ '4402 dS 

-■3478 dt 


III. Occultations during 

the Lunar Eclipse, 

1899 December 16. 

1. 

- 3-24 = ds 

-•4330 dp 

4- '6954 da 

4 " '6501 dd 

+ -2544 dt 

2. 

— 277 = ds 

-•4950 dp 

4- *5822 da 

+ 7711 dd 

+ -1988 dt 

3- 

- 3‘59 = ds 

— *2908 c/p 

4- *8910 da 

4- *2406 dd 

+ '343 1 dt 

4- 

+ 0-26 = ds 

— •0643 dp 

— 7082 da 

4 - *4364 dS 

-•4354 dt 

5- 

+ 2-52 = ds 

— •1817 c/p 

— •4068 da 

4- *8982 dd 

— -2492 dt 


It was obvious, of course, that the above reductions were 
altogether insufficient in number to give by themselves any 
reliable value of the semi-diameter of the Moon, but an attempt 
was made to examine in a general manner the results which 

N 2 
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156 Mr. Plummer, Oecultations of Stars etc. lxi. 3, 

might be deduced from the equations. The attempt was not 
successful, and it is not without interest to consider why this 
was so. The case of the ten disappearances of stars in the 
Pleiades, I. 1-10, is typical. The equations were given equal 
weights after division by the coefficients of dt, which term was 
afterwards neglected of course. The star places are known to be 
very accurate, and the residual, therefore, is fairly attributable 
to the error in the observed quantity—the time, and can be 
treated in accordance with the theory of least squares. Theoreti¬ 
cally, this method seems reasonable, but it may have the 
result of giving an overwhelming weight to the result of an 
occultation at a point of elevation or depression in the Moon’s 
figure. It would seem, then, that a compromise is desirable in 
some cases, but the question is one which I do not pretend to 
decide. As regards the equations referred to, this difficulty did 
not exist, and the four normal equations were at once formed in 
ds, dp, da and d8. It was then found that these equations were 
not, in any real sense, consistent and independent, as would be 
necessary for the determination of the four unknowns. In fact, 
they could be reduced to two without sacrificing any of the infor¬ 
mation which they were capable of giving. This result can be 
referred, I think, to two simple causes. In the first place, the 
observations do not extend over a long interval of time, and the 
Moon’s change of position is small. Now in any position it is 
clear that if a compensating change were made in the radius of 
the Moon, the Moon’s distance from the observer might be 
changed without any effect on the phenomena. Consequently 
the Moon’s apparent semi-diameter and parallax cannot be de¬ 
termined independently from oecultations observed during a 
short time at a single observatory. In a similar way if dis¬ 
appearances alone are considered, it is evident that the effect of 
moving the centre of the Moon along its own apparent path 
while changing at the same time its semi-diameter so that the 
advancing point of the edge retains its position, will be very 
small in the case of occulted stars which lie near that path. In 
fact, in both cases we shall be trying to determine small quanti¬ 
ties of the second order, and it is inevitable that the results will 
be expressed in terms of the residual errors. The difficulty can 
only be overcome by combining results obtained at distant ob¬ 
servatories, and it is in the hope that the reductions given above 
may be utilised in this manner that this work has been under¬ 
taken and published. 

Oxford : 1901 January 7. 
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